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We study physical properties of a Luttinger liquid in a superlattice which is characterized by 
alternating two tunneling parameters. Employing the Bosonization approach, we describe the corre- 
sponding Hubbard model by the equivalent Tomonoga-Luttinger model. We analyze the spin-charge 
separation and transport property as the difference between the two tunneling parameter increases. 
We suggest that cold Fermi gases trapped in a bichromatic optical lattice and coupled quantum 
dots offer the opportunity to measure these effects in a convenient manner. We also study the clas- 

C/3 , sical Ising chain with two tunneling parameters. We found that the classical two-point correlator 

uJ ■ decreases as the difference between the two tunneling parameter increases. 
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I. INTRODUCTION 



Quantum many-body systems of one-dimensional interacting fermions have attracted enormous interest for more 
than 5 decades. Contrary to what happens in two and three dimensions, these systems cannot be described by the 
Landau theory of normal Fermi liquids. The appropriate paradigm for ID interacting fermions is instead provided by 
the Luttinger liquid concept introduced by Haldane. The distinctive feature of the Luttinger liquid is that its low- 
energy excitations are collective oscillations of the charge or the spin density, as opposed to individual quasiparticles 
that carry both charge and spin. This leads immediately to the phenomenon of spin-charge separation, i.e the fact 
f^ that the low energy spin and charge excitations of ID interacting fermions are completely decoupled and propagate 
psj I with different velocities. Despite the firm theoretical basis upon which the Luttinger liquid theory rests, there has 
C^ been precious little compelling experimental evidence that real one-dimensional electron gases are anything but Fermi 
^/-^ , liquids. In recent years, it has become possible to fabricate single channel quantum wires but unwanted impurity 
f*^ ' causes backscattering and localization, thus destroying the Luttinger liquid phase. Fortunately, there is another 
f^ , experimental system which is expected to exhibit Luttinger liquid behaviour, and does not suffer from complications 
^~~^ • associated with impurities-namely one-dimensional quantum Fermi gases in optical lattices. 

^ [ In this context, it has now become possible to trap ultra-cold quantum gases in quasi-lD optical lattices. Much 

■ ^ of the theoretical work has been on ID Bose gases but more recently Recati et al [l| have studied one-dimensional 

^% quantum gases of fermionic atoms in optical lattice using Luttinger liquid approach while Polini et al. Q] have 

studied spin-drag and spin charge separation of atomic Fermi-Dirac gas in a one-dimensional optical lattice. Using 
superposition of optical lattices with different periods [Sj , it is now possible to generate more sophisticated periodic 
potentials characterized by a richer spatial modulation, the so-called optical superlattices. An important and exciting 
application of optical superlattice is quantum computation 4] . The physics of one-dimensional optical superlattices 
has been a subject of recent research, including fractional filling Mott insulator (MI) domains [5], dark [6.1 and gap 
[31 solitons, the Mott-Peierls transiti on pi , non-mean field effects H, phase-diagram in two colour superlattices [10| . 
Bloch-Zener and dipole oscillations [11|, collective oscillations [l2] and Bloch and Bogoluibov spectrum [1^. In 
section II, we discuss the influence of such a kind of superlattice structure on the spin-charge separation and the 
compressibility of the Fermi gas. 

Often referred to as artificial atoms, semiconductor quantum dots offer an unprecedented possibility of constructing 
at will and exploring situations ranging from practically single atom to a fully solid state many-body systems. The 
nanofabrication possibilities of tailoring structures to desired geometries and specifications, and controlling the number 
and mobility of electrons confined within a region of space, makes these structures unique tools to study transport 
properties. Quantum transport in arrays of tunnel coupled quantum dots have attracted attention for the past few 
years [14l4l8l |. The controllable quantum properties of the electron in such systems opens the possibility of their 
application to schemes of quantum computers !19;]. In section III, We study the conductivity of an electron that 
experiences an asymmetric tunneling when going to the left and right in an one dimensional tunnel coupled quantum 
dots. 

The spin 1/2 Ising chain is considered the prototypical system for quantum phase transitions 12^1. However very 



development of molecular beam epitaxy, it is now possible to envisage a superlattice in which the exchange constant 
varies from layer to layer. Very often one finds interesting properties in these systems [2l|. Magnetic excitations in 
superlattices were studied in numerous works (see [22] for a brief review) . Hinchey and Mills [23] have investigated a 
superlattice structure with alternating ferromagnetic and antiferromagnetic layers. A common feature that connects 
the three systems studied in this work is the asymmetric tunneling. 

II. THE HUBBARD HAMILTONIAN AND THE EQUIVALENT TOMONAGA-LUTTINGER MODEL 

Our aim will be to study strongly correlated systems in one spatial dimension. These are typically systems of 
interacting electrons but we will be interested in cold Fermi gases also. The prototypical interacting electron system 
is the Hubbard model. This is the lattice model whose Hamiltonian in one dimension is 



Hfh = - V J," (c]+i_^cj,^ + 4,aCj+i,a) +uY^ ^i,t"j,; (1) 
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The first term describes the hopping process, in which an electron can move from one site to the next with site 
dependent amplitudes Jj which takes two distinct values, Jq + (— l)-* — ■ The hopping process preserves the spin 

projection a. For cold Fermi gases, such a kind of hopping terms can be created by superposition of two optical lattice 
of different periodicity (as described in section 3) while for electron system, coupled quantum wells with appropriate 
voltages can generate such a hopping term (as described in section 4). The second term describes the local Coulomb 
repulsion (C/ > 0) between opposite spin electrons residing on the same site. For cold Fermi gases this would be the 
two body interaction as discussed in the next section. The Cj^o- operators are the usual annihilation operators with 
anti-commutation relations. Also, hj^^ — cj „Cj_„ is the Fermionic number operator. In the next section, we will show 
that the Hamiltonian for the cold Fermi gas in an one dimensional bichromatic optical lattice can be reduced to the 
Fermi-Hubbard Hamiltonian [1] In this section, we will derive the equivalent bosonized Hamiltonian of the Fermi- 
Hubbard Hamiltonian. The technique of bosonization is a powerful tool to study the spectrum of low-lying excitations 
and correlation functions of one-dimensional systems. Let us first look at the non-interacting limit ([/ = 0). In this 
case, the Hamiltonian can be easily diagonalized by means of Fourier transformation. We define 

^^.- = ;^ E {'lo - ^(-1)^4..} e^'"'^' (2) 

Here L is the number of lattice sites with periodic boundary conditions Cj+^.cr = ^j,a- As the fermions move from 
one well to the next, it acquires as additional phase, which depends on the height of the barrier. As the height 
alternates, the phase also alternates. This picture is conveniently represented by the j dependent inverse Fourier 
transform of equation (2). Substituting equation (2) in the non- interacting fcrmi-Hubbard Hamiltonian, we have 

H = -2jocos2kd{4lcl^ - cllcl^} - Aosin2fcd {c^,|,4,, + cllcQ (3) 

Finally, H can be brought to diagonal form by defining operators 
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1 Jocos2kd 

Jk,a = a;k,a (4) 



Cfc 



1 Jo COS 2kd 

O fk.a = I3fk,a (5) 



Where \ak\^ + |;3fcP = 1 and Ck = J W^ cos^ 2kd + Ag sin^ 2kd. This yields: 



The ground state for N fermions corresponds to filling up all the states, from the lowest energy up, until the A^ 
lowest-energy orbitals are filled up (taking into account spin degeneracy). The highest occupied level is the Fermi 
level, its energy the Fermi energy Ep and its wave- vector the Fermi wave- vector kp. The relationship between N 
and kp is N ~ 2kpL/TT or n = N/L = 2kp/TT. When we take into account interactions and if U << Jj (perturbative 
region), it is natural to assume that only low energy states will be affected. This is reasonable within second order 
perturbation theory. We now introduce in the usual way the right movers (around +kp) and left movers (around 
— kp). We then have two linearized spectrum around the two Fermi points 

e(fc) ~ vp{k — kp) =^ Rightmovingbranch (7) 



e{k) — —vp{k + kp) ^=> Leftmovingbranch (8) 

according to the sign of the velocities. The relationship between this spectrum and the lattice one is given by: 



^F = ^, \k=kF = — . = yd) 



de{k) _ rfsin(4fcFd)(4Jo^- Ag^ 

cii, \k=kF I 

"'^ h^(4:J^cos^2kFd + Alsm^2kpd) 

Note that, we effectively restrict ourselves to low energies. Corresponding to the right and left moving fermions, 
we can introduce fermion annihilation(creation) operators cl„, where v = R,L and the respective density fluctuation 

operators Pi^.a{q) = ^k'^t ai^ ~^ l)'^i';<^i^)- Note that p,y^cr{—q) — ptail)- These particle-hole excitations can be 
written in terms of bosonic creation and annihilation operators: 
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The normally ordered number operator is defined as: 

N,,a = Y. ■ 4,a(fc)c.,.(fc) ■.,iy = L,R (11) 



and 



^-A=c,s = -^(iV,,t±A^.,i) (12) 



As a part of the bosonization process, we also introduce boson field operators: 



R,L,a - ^f Z^ _e-"9/2 [e^ii-bf^ - e^^'i'-bf^^) (13) 

^^ q>Q V9 



(l),..\=c,s = -y= {(f>i^.ti^) ± (f>i^.li^)) (14) 



CA=.., = 4(Ct±c,.) (15) 



H^Ho + Hm (16) 
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B^r^t - ^E {1 [^ (^^A + ^Ia) + 92X^rxNl^ ]+^T.f^\ ^ I ^ E ^ (^-'^-a)' : -92X : {d.<f>Rx) {B^lx) : 

(18) 
Where, gi.x=c.s ~ ffiji i5i,±7 i = 2,4. g^ ii is the intraspecies interaction which is zero for cold Fermionic atoms. 
gi^± — Ud is the interspecies interaction for cold Fermionic atoms. 52A is the strength of the forward scattering 
between particles belonging to different branches and with different spin state, while g^x is the strength of the forward 
scattering between particles belonging to same branches and with different spin state. We now have two decoupled 
sectors corresponding to charge and spin excitations. The Hamiltonian in equation (), can be diagonalized by the 
following Bogoliubov transformation: 

dl^ = cosh 7S«, + sinh 75^^3, (19) 

rf2, =sinh76^;, + cosh7S^^;, (20) 

Q2X 

Where, tanh27 = A = . Finally, the diagonalized Hamiltonian is: 

2ttvf + gix 

H = hJ2uxY. I^gld'^x + JZ (''"^^' + ''^^"^') ' ^ = 1' 2 (21) 



Where, 



Nx^Nrx + Nlx, (22) 

Jx - Nrx - Nlx (23) 



Ux 

VNX = — 24 



vjx ^ uxgx (25) 



= / 2^^F+g4A-g2A (26) 

^' 27:vF + gix +92X 



In the above discussions we have neglected the back scattering and Umklapp terms. For the case of cold Fermi 




Figure 1: Schematic drawing of the optical superlattice with alternating big and small wells. 
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Figure 2: Plot of normalized compressibility as a function of a for U/2Jo = 0.2 and Kpd — 0.1 



III. THE HUBBARD MODEL FOR A COLD FERMI GAS IN A BICHROMATIC OPTICAL LATTICE 



We consider an elongated cigar shaped (quasi-lD) dilute gas of fermionic atoms of mass m with two internal ground 
state, \a —'t,l>, representing a spin-1/2 system. We will assume that the two internal levels are equally populated i.e 
N^ = N^. The atoms are cooled below the Fermi-degeneracy temperature ksTp ^ Nhuii, N = Nf + N^ is the total 
number of particles. The condition for a quasi-lD system is a tight transverse harmonic trapping with frequency uj± 
exceeding the characteristic energy scale of the longitudinal motion. In this way the transverse degrees of freedom 
are frozen. Because of quantum degenaracy, the longitudinal motion has all the energy levels up to the Fermi energy 
ep ~ ksTp filled. Typical values of u!± and uii (frequency of longitudinal confinement) are in the range 27r(300 — 400) 
Hz and 27r(2— 10) hz respectively. Thus we require the total number of particles to be restricted by A^ < ujj_/uji, which 
is typically of the order of few hundred. Because of the Pauli principle, at low temperature, only s-wave collisions 
between atoms in different internal states are allowed. Therefore, all the relevant interactions are characterized by 



inter component scattering length a. The strength of the effective ID interaction is g — 

is the harmonic oscillator transverse length and a < l±. 

Thus the system is described by the following one-dimensional Hamiltonian 



27rh^ 



-, where l± — 



muj± 



H = ^ dx^jlix) (-^dl + V,,t{x)\ ij, + g I dx^\{x)^l{x)^^{x)^^{x 



(28) 



Where ipa is the ID field operators for atoms in state a. The external potential Vext{x) — Vl{x) + Vop{x), 
includes both the longitudinal confinement Vl{x) = —mujfx^ and the two-colour optical lattice potential Vop{x) = 
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Figure 3: Plot of the charge velocity u^ (thin curve) and spin velocity u^ (thick curve) as a function of kpd (left plot) for 
a = 0.1 and U/2Jo — 0.2 and as a function of a (right plot) for kpd — 0.1. 

particular case d2 = 2di — 2d. In addition, we will consider the case when the optical lattice dominates over the 
harmonic potential. We expand the atomic field operators in the lowest-band Wannier basis 



tpcr — 2. W{x — Xj)c 



Jif^' 



(29) 



Where, W{x — Xj) is the Wannier function centered at the j'*'' site and Cj^a- is the annihilation operator for a fermion 
in the j*'' site with spin cr. Sunstituting equation (2) into equation (1) and retaining only the nearest neighbour terms, 
we get the Fermi-Hubbard Hamiltonian: 



H 



FH 
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(30) 



Here, the onsite energies are taken to be zero. Jj is the site dependent tunneling and takes two distinct values, 

Jo + (— l)"* — . The strength of the effective on-site interaction energy is U = g J dx\W (x)]"*" . 

The simple form of the Hamiltonian [2T] makes the calculation of some physical properties rather straightforward. 
One important quantity, the compressibility of the gas is written as: 



K 
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(31) 



Here, kq is the compressibility of the noninteracting gas. The density of the gas is approximated as homogeneous. 
This is true if the trapping potential along the optical lattice axis is very weak and the also the depth of the optical 
lattice is small. The systems interacting with replusive interactions, the optical lattice reduces the compressibility of 
the system, since the effect of repulsion is enhanced by the squeezing of the condensate wave function in each well. 
Uc/gc fixes the energy needed to change the particle density. From figure 2, we note that the compressibility decreases 
with the parameter a. This means that with increasing strength of the secondary lattice, the energy required to 
change the particle density increases. Figure 3 shows a plot of the spin velocity (w^) and charge velocity (w^) as a 
function of the quasimomentum (kpd) and the parameter a. As the condensate moves across the Brillioun zone, the 
difference between the spin and charge velocity increases and is maximum at 7r/4. This suggests an effective mechanism 
to observe the spin-charge separation. One can move the condensate across the Brillouin zone by accelerating the 
condensate in the optical lattice. On the other hand the spin-charge separation decreases with increasing a. This is 
perhaps due to the fact that as the strength of the secondary lattice increases, the condensate becomes more localized. 



IV. LINEAR ARRAY OF TUNNEL COUPLED QUANTUM DOTS 



r~^ r~> r~> r~> r~^ r~^ r^ r~> 
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Figure 4: Schematic drawing of the chain of tunnel-coupled quantum dots. Notice the alternating tunneling coefficients. 




Figure 5: Normalized conductance versus a for U/2Jo — 0.2. 



(GaAs/AlGaAs). The individual tunneling rates are determined by the voltages applied to the gates defining the 
corresponding intcrdot tunneling barriers. Thus appropriate gate voltages have to be applied to reproduce the 
Hamiltonian [1] 

The conductance G is found as 



G = Go/v/l + 7 



(32) 



Where, Gq — 2e^//i and 7 — 2U/{Trhvp). In figure 3, we plot the dimensionless conductivity {G/Gq) for repulsive 
interaction versus the parameter a = Aq/2Jq. We notice that as a increases, the conductivity decreases which again 
is an indication of localization of the electrons. The correlation functions give information about the tendency of the 
system to show long range order. The correlations that decay the slowest are the dominant ones. The bosonization 
method makes the calculation of the correlation functions straight forward. Here we will focus on three most dominant 
correlations, i.e., oscillatory part of the spin density (spin density wave, DJ^^, ~ l/x^'''="'"^=\ ^^dw ~ l/x^^"'^^^^"^) and 
oscillatory part of the charge (charge density wave, Dcdw ^ l/x^-^"'^^''^). An analysis of the factors gc and gs reveals 
that as a increases the only correlation that becomes dominant is D^^^, i.e the transverse component of the spin 
density wave. 

V. THE CLASSICAL ISING CHAIN 

It is known that the quantum transitions in the quantum Ising model in d dimension is intimately connected to 
finite temperature phase transitions in classical Ising model in D — d + I dimension. The D = 1 and A^ = 1 classical 
Ising model does not show any phase transition but it has regions where the correlation length becomes very large 
and the properties of these regions are very similar to those in the vicinity of the phase transition points in higher 
dimensions. Here we will consider the D — 1 and A'^ = 1 classical spin ferromagnet, more commonly known as the 
ferromagnetic Ising chain. This chain has the partition function 



^= X! ^^Pi^H), 



(33) 
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Figure 6: Two point correlator as a function of lattice site j for three values of Ao/2Jo = 0.01, 0.25, 0.5 
. The thickness of the plots increases as the value of Ao/2Jo increases. 
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Here M is the total number of Ising spins and h is the external magnetic field. We will assume periodic boundary 
conditions, therefore cr^j,^ = a^ . Now following Ising, we write Z as a, trace over a matrix product with one matrix 
for every site. 



M 
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where 



Ti,i{ai,a. 



i+l) 
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and 
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Ti^i will have values different for i even and i odd. The matrix T1T2 is identified as the transfer matrix of the Ising 



chain. Let us now define a matrix T3 = Ti.i^euenTi i=odd- ^^ ^^^ limit il/ — > 00 one can show that 

2^2^ iM/2 _ A//2 A//2 



Z = Tr[TlT^Tz] 



+ e. 



(38) 



where ei and €2 are the eigenvalues of T1T2T3. For the case h — (no magnetic field), T2 — 1. The eigenvalues are 
found as 



2cos/i(Jn) , , . ^ , —-- — 7 

e-i = , — ; ' , , ^sinh{2Jo) + sinh{Aa) 

y/ sinh[2Jo + Aq) 



(39) 



The eigenvalue ei increases with increasing Ao and saturates at 1.2 at Sq — 2Jq while the other eigenvalue £2 goes 
to zero at Aq = 2Jo. 

Now we calculate the correlation function exactly. For simplicity we consider the case of zero external field {h = 0) 
and describe the two-point correlator 

<<'^I>=^E^^"''^<'^|- (41) 

In the limit of an infinite chain, the two-point correlator in terms of continuous variables is derived as 

C{t, 0) =< cr(T)a(0) > e-l^l/«, (42) 

where the correlation length ^ is written as, 



-z — — { ln[coth{Ja)] -\ — In 
t. a I 2 



sinh{2J()) 4- sm/i(Ao) 



sinh{2Jo) — sinh{Ao) 



(43) 



Here, r = ja, a is the lattice spacing. Figure 6. displays the two point correlator as a function of lattice site j 
for three values of Ao/2Jo = 0.01,0.25,0.5. Clearly, we see that as Ao/2Jo increases the correlation decays faster 
indicating the fact that the spins are getting localized. 

VI. CONCLUSIONS 

In conclusion we have studied three one dimensional superlattice systems (characterized by two tunneling param- 
eters) namely, atomic gases in one dimensional superlattice, linear one dimensional array of quantum dots and the 
one dimensional classical Ising chain. In particular for atomic gases, we found that as the difference between the two 
tunneling parameter increases, the difference between the spin and charge velocities decreases. This is attributed to 
the increasing localization of the atoms in the wells of the optical lattice. On the other hand for the case of linear array 
of quantum dots, as the difference of the two tunneling parameters increases the conductance decreases attributed to 
the pinning of the electrons. For the classical Ising chain, the two-point correlator decreases with increasing strength 
of the superlattice which is attributed to the localization of the spins. In general we conclude that atoms, electrons 
and spins are comparatively more localized in a superlattice structure. This study demonstrates that by tuning the 
two tunneling parameters, one can coherently control the transport properties of a superlattice structure. 
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